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I.  Abstract 


Vh  =  {«  e  W  :  „({,,)  =  E"  «*,(«)*,(.,)} 


N\  AR  3  0  A990 

l 


(3.13) 


A  feed-back  extension  procedure  is  developed  for  the  numerical  solution 
of  a  clast  of  nonlinear  boundary  value  problems  associated  with  anti-plane 
shear  or  Hencky's  theory  of  plasticity.  This  extends  previous  results  using 
dimensional  reduction  tu  energy-asymptotic  format.*”^  (  j 

2.  Introduction 


In  an  earlier  paper  |4|,  the  method  of  dimensional  reduction  for  quasilinear 
boundary  value  problems  was  introduced.  A  generalisation  was  proposed 
which  allows  for  the  possibility  of  different  order  of  dimensionally  reduced 
models  in  different  parts  of  the  underlying  domain.  This  paper  is  an  at¬ 
tempt  to  fulfill  that  promise  with  the  purpose  of  making  the  method  of 
dimensional  reduction  still  more  efficient  and  robust. 

As  in  (4|,  the  basic  idea  is  to  find  a  minimiier  uN  of  the  given  energy 
functional  in  a  proper  subspace  VN  which  is  characterised  by  the  basis 
'unctions  {^,){t„: 

Vh  =  {E7,»c,|f)*,(«l)} 

where  £  =  G  |0,  1 1,  r;  =  xj/d,  €  |-d,dj,  and  d  denotes  the  half  thick¬ 
ness  of  the  domain.  Thus  the  model  of  order  N  of  reduced  dimension  was 
introduced.  See  [4|  for  the  choice  of  {0,}yLo  a^d  related  convergence  prop¬ 
erties  (optimal  rates)  as  d  1  0  or  N  —  oo. 

Due  to  the  singularities  which  can  stem  from  the  loading  or  the  presence 
of  corners,  it  is  necessary  (for  efficiency  and  accuracy)  to  be  able  to  intro¬ 
duce  higher  order  models  near  these  layers  only.  In  this  paper  we  propose 
a  feed  back  extension  procedure  that  facilitates  this  by  allowing  different 
orders  JV,  in  different  parts  of  |0,l). 


3.  Notation  and  Model  Problem 


We  shall  confine  our  study  to  the  following  class  of  problems.  Find 
such  that 

Vu  €  W ,  Au(v)  =  G(v) 

u  €  IV 

(3.1) 

where  r 

Au(v)  =  /  F(|  V.jU  j2)V,|U  V.iv  d£drj 

Ml 

(3.2) 

G(v)  =  d1-/*  0K)(vtf,  1)  +  vtf,  -1)|  dt 

4 n 

(3.3) 

=  l  +  r,  n  €  N,  t  €  K  \  R, 

(34) 

[  €  R  characterises  three  asymptotic  ranges  of  loads  : 
p  <  -  such  that  the  limit  traction  on  T±  as  d  J  0 

(  is  0  for  ix  <  0,  finite  for  p  =  0,  and  infinite  for  p  >  0 

(3.5) 

n  =|0,t|x|-  1,1| 

(3.6) 

r«  =  ((o)x|-i,i|)u((i)x|-i,  ■() 

(3.7) 

rt  =|o,i|x  (i) 

(3.8) 

r.  -|o,i|x  (-i) 

(3.9) 

W  -  W|‘v-3"+I(n)  =  (u  e  vv,  a~+3(n);  v|r„  =  0) 

(3.10) 

«  ,  a  >  a  , 

■*“  i* 

(3.11) 

This  scalar  problem  corresponds  to  finding  a  minimiier  in  W  for  the 
energy  in  anti-plane  shear  in  finite  elasticity,  |6|  and  |7|  and  the  torsion 
problem  for  a  bar,  see  |8|  and  |9j.  See  |4j. 

We  define  the  dimensionally  reduced  solution  of  order  N  to  be  the  so¬ 
lution  us  in  Vn  C  W  for  which 

VvCVnCW,  Au*(v)*C(v)  (3.12) 

given  Vn  a  subspace  in  W  of  the  form 


'Partially  ■upporttd  by  the  OlHce  of  Naval  Research  Contracts  No.  NOOH-tS-K-OIGO 
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The  family  of  subspaces  (Kv)jv,,u  “  characterised  by  the  choice  of 
called  the  basis  or  Ansats  functions. 

In  (4|  these  basis  functions  were  selected  to  yield  the  optimal  rate  of 
convergence  of  |(u  -  un((h>  14  d  ;  0.  We  thus  had  to  select  0y  to  be 
a  polynomial  of  degree  2j.  Importantly,  the  tame  choice  is  valid  for  all 
three  ranges  of  loads  (three  signs  of  n  in  (3.5)  and  (3.3)).  For  F  in  (3.4) 
depending  on  rj,  (4|  indicated  that  the  same  procedure  would  yield  0y  to  be 
a  nonpolynomial  solution  of  a  second  order  Sturm  Liouville  problem.  (See 
also  |5,  Remark  3.9|.) 

Let  un  be  the  Nth  partial  sum  in  the  formal  asymptotic  expansion 
as  given  in  |4|.  Let  Di  be  the  operator  defined  by  Dju  =  ^  mapping 
L>om(Vt)  =  W,'3"+5(0,  l)nr0l  h”  -  £»»+»( o,  1).  We  cot  for  p  <  0. 
Theorem  S.l  Let  jj  <  0  and  n  e  Z.  Letups  W‘“  ke  hounded  there 
independently  of  d.  Let  0  6  Dom(D **)•  Then  there  ezijte  Cft  independent 
of  d  tuck  that 

||u-u„ll„, 

Since,  for  a  given  practical  problem,  we  cannot  depend  on  d  being  suf¬ 
ficiently  small  to  ensure  that  a  given  tolerance  criterion  can  be  satisfied 
via  the  previous  theorem,  we  have  considered  in  |4j  to  increase  N.  Again, 
optimal  rates  in  this  scenario  (d  fixed,  N  increasing)  were  established  in  (4[. 
FYom  the  computational  experience  in  |4|  and  elsewhere,  it  became  clear 
that  it  was  unnecessary  (read:  wasteful)  to  increase  N  uniformly  every¬ 
where  in  (0,l|.  Rather,  there  were  clearly  defined  layers  (near  the  boundary 
and/or  rough  spots  in  the  load).  We  propose  to  increase  N  near  these  layers 
only  as  our  extension  procedure. 

Lei  I  =■  (0.1)  =  lC,I„  end  L  n  I,  =  |,i  /  ,  Vi,;  e  (l,m|.  Let 
N  *  (^.)"  I  be  an  m-vcctor  of  nonnegntive  integer*  (Af,  =  no.  of  baeie 
function!  ueed  in  It).  Conaider 

vn  -  {v  :  »(?,!»)  “  ^J'.0vj(0^y(i|)  auch  that  (3.14) 

*  =  PIU.  v»(‘)  *  0  for  f  e  U;>/v,I,} 
a  subs  pact  of  VN.  Solving 

V»6V*  cj»,  .<«*(*)  wC(u)  (3.15) 

for  v  €  V x  i*  the  generalised  dimensionally  reduced  Calerkin  problem. 

A  key  ingredient  in  the  selection  of  the  distribution  of  orders  M  -  the 
local  a  posteriori  estimators  -  will  be  developed  in  the  following  section. 

4.  Local  A  Posteriori  Error  Estimators 

Define  the  eelimatorfor  (0,1)  and  order  M  as 

*  I  j^lU’ini  (<•!) 

where  e  €  fl )  is  the  solution  of 

Vv€ff‘< 0,(0):  f_J  l£ij^^  =  CH-/tun(n)  (4.2) 

the  right  hand  side  being  the  residual  (Au-  Au*  )(v).  Although  e  is  not  well 
defined,  £*  and  Btt(M)  are,  provided  the  following  solvability  condition  is 
satiafied 

Vc  €  H'{ 0,  l)  :  ]  3i‘)2c[()d'-»  =  (4.3) 

f'  J  F{ \VJV,  l)’)a"~c '(flddfdo 

'-I^O 

However,  this  is  satisfied  (even  for  c  €  IVp**"4’*)  if 

lespon((^),NMO)  (4.4) 

cf.  (3.12)  and  (3.13).  This  condition  is  met  for  any  choice  of  basis  functions 

with  optimal  rates,  see  |4|. 
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Similarly  define  the  local  error  estimator 


£*«<(*)  =  I ['t  (~)a  ddl-dr,]*,  I  S  i  <  m 


«•*) 


Lemma  6.2  There  exist*  a  constant  matrix  A  such  that 


At  in  |2|  we  define  upper  (lower)  error  eetiniator  to  mean 
||«  -  Wjw  ||w  <  (>)£st 

Theorem  4.1  Let  a  and  u*  he  the  exact  and  dimensionally  reduced  so¬ 
lutions  (see  (3.1)  and  (3.15)).  Then  Est  as  defined  in  (l-l)  if  an  upper 
estimator,  i.e. 

l|u  - 

Proof;  Bound  from  above  and  below  (Au  -  Au*)(u  -  u* ).  j  j  j  j 

In  the  language  of  (3|,  Eat  ie  a  guaranteed  U-estimator  (0-«»timator). 
Another  attractive  property  of  Est  it  that  £jf/||u  -  tends  to  1  as 

d  tends  to  sero  for  ft  <  0  and  fl  sufficiently  smooth. 


Theorem  4.2  Let  u  and  u m  he  the  exact  and  dimensionally  reduced  solu¬ 
tions  (cf.  (3.1)  and  (3.15))  and  have  gradients  hounded  uniformly  in  d.  Lei 
/4  <  0  and  fi  €  Dom{D {*).  Then  Est  as  defined  *n  (d-1)  i*  asymptotically 
exact: 

£ft(/f )  =*  ]|u  -  ||h»(I  +  0(d)) 

Proof:  Due  to  restrictions  in  length,  we  merely  mention  that  one  can 
establish 

<  |«  -  u*  fl# .  -M^-(i  +  0(d)) 

and  bound  the  middle  factor  on  tlii.  right  hand  aid.  from  above.  .... 


6.  Computational  Agpectg 


where  </r+i,  cn+3  were  defined  in  the  prertoue  Lemma. 

Proof:  Consider 

~LU“  =  liph  =  +  «!*.*>«♦>) 

We  denote  by  (■,  •)  the  inner  product  in  L3(0, 1).  Now  let  v  *  4>0 1,  then 
0  =  (/Sd1-**.  (^o(l)  +  ^o(-l)*)  -  1^0*) 

=  (2/M'-«  .)  +  ({Lu„ .  do),  a)  -  (/•(|VduJ,  |J)  ™  |L,.  *) 

=  (2^d1-'‘-8llew  +  1-fl1,«N+J-P(|t7J„*|’)a-u*|L1,*) 

dtj 

Next  letting  v  =  fix,  we  obtain  similarly 

0  =*  — (Bile***  +  +  F(|Vjujf  x) 


where 


Bn 


r 


iil 

<P  dri 3 


ts+sti-  t  d  di) 


for  1  <^jj .<_?•  The  matrix  B  is  invertible  into  A. 


fill 


We  next  introduce  the  heuristic  principle  which  will  guide  us  to  an  efficient 
extension  procedure  based  on  the  local  a  posteriori  error  estimators. 


Heuristic  6.1  Let  the  error  associated  with  the  generalised  dimensional 
reduction  he  estimated  hy 


FYom  an  implementational  point  of  view,  the  nice  mathematical  properties 
of  Est  established  in  the  previous  section  will  not  suffice,  since  findiug  e  as 
a  solution  of  a  second  order  O.D.E.  might  be  too  costly.  Therefore  we  give 
now  some  formulae  that  can  be  used  to  compute  e  and  Est  in  practice. 

First,  let  us  introduce  a  basis  in  La(-1,  1)  with  which  we  will  work: 


and  the  cost  (work)  he  estimated  hy 

ERi*  (*<•*) 


do(n)  *  1  (5.1) 

W*l)  ■  /  <J-|(0  (it  for  ;  >  1 

-i 

where  1,  is  the  jth  Legendre  polynomial. 

Lemma  6.1  The  following  formulae  for  e  defined  in  fflj  hold: 

,  =  E*+Va 


Then  we  atm  al  achieving 

EstUNi)  -  Est^Si  -  1)  «  V (M,!,)  -W[Si  -  l.L) 

hy  increasing  W,  hy  1  where  the  error- cost  quotient  is  maxtmai. 

Reasoning:  Minimising  the  error  at  fixed  cost  with  respect  to  Nit  yields 
via  Lagrauge's  multiplier  and  a  backwards  difference  approximation  the 
proportionality  aimed  at  in  the  Heuristic.  , 


and  furthermore 


A  typical  choice  for  workestimate  is 


0  for  1  <  j  <  N 

<?  It ,  CuM  d,  dn/  /!,(dj)a  d<,far,-N  +  l.fl  +  7 


where  Lis  ■» 


Proof;  e  €  V/#+a  so  the  first  assertion  is  dear  and  Jj  «  Ej^+,3^'  (fj)ey  and 
if  we  denote  by  (,)  the  inner  product  in  L3(- 1, 1)  we  have  V*  € 


ljTn’*>)dn 
I  y)*i  >*)*)  *1 

d(Mu-^u*l(d,*) 

0  for  1  s  }  <  N 


Hence,  e,  ■=  0  for  1  <  j  <  N  and  the  two  remaining  components  are 
obtained  in  a  similar  fashion  using  now  the  definiton  of  A  in  the  next  to 
last  line  and  the  fact  ^,(l)  «■  ^,(-1)  «■  0  for  j  >  1.  mj 

The  following  formulae  derived  from  the  previous  ones  art  still  more 
useful: 


^(^.L)-(oi^  +  l)-|L|  (5.2) 

for  some  choice  of  positive  a,  ,  in  1,2. 

Note  that  we  selected  the  prime  functions  of  Legendre  polynomials  as 
basis  merely  to  be  able  to  establish  computational  formulae;  a  change  of 
basis  within  the  same  span  merely  requires  a  linear  transformation  in  order 
to  modify  the  formulae  for  the  new  choice  of  basis  functions. 

FYom  the  computational  point  of  view  it  is  rather  important  exactly 
which  basis  functions  one  selects  (this  has  to  be  done  hierarchically). 

8-  Choice  of  BaeU  Function* 

Ut  -  (dj)f.o  »»d  v  -  («,){'.„  inch  that 

=  U 

Th«  f entrained  G.lerkin  problem  (3.15)  tr. inform,  to  tk.  following  .y.t.m 
of  t>  D.  E.e: 

-  4-d{P{u,  v')V)  +  i<j(y,  u')u  -  R  (e.i) 

dx  a 

where  P  and  Q  are  matricee  defined  by: 
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Plj  -  /‘  /•(|VJu»|s)*0J  4|  («.2) 

■'-1 

<3.7  -  /‘i/,(|VJu„|>)lo;*l  (6.J) 

for  1  <  I,;  <  TV.  Sine*  the  system  (6.1)  is  hard  to  analyse  in  its  nonlinear 
form,  we  will  bracket  with  linear  ones.  If  l|Vju*|Joo  ^  Af,  then 

C/')  <  <  (1  +  M2n)[PtinV'lU') 

(Q,tnU\  U‘)  <  (Q(£/,t/')£/\£/')  <  (1  +  M2n)[QliHUltUt) 

where  /*'"  and  QUn  are  defined  as  in  (6.2)  except  with  F  si.  This  allows 
us  to  analyse  some  of  the  behavior  of  (6.1). 

An  elementary  Saint- Venant  like  principle  holds  for  a  related  linear 
boundary  value  problem  posed  over  the  semi-infinite  strip:  f)°°  =]0,  oojx  ]  - 
1,  l|  with  boundaries  Tg8  *  {0} x)  -  1,  l|,  and  TJ*  which  is  defined  analo¬ 
gously  to  (3.8)  and  (3.9)  respectively.  The  function 

i  .  „oo  fcir* 

u(x,  y)  =  Efc-o ak  cos  —  exp  -  — 

is  the  solution  of 


A,iu  =*  0  in  n®°; 
du 

—  —  0  on  T?; 
du 

u  =  g  on  rg° 

Here  A*  =  k  €  No  are  the  eigenvalues  corresponding  to  the 

eigenfunctions  given  by  the  following  B.V.P.  (in  the  y-directiou) 

**  0  in  )  -  d,  d| 

*  0  at  ±  d 

The  eigenvalues  may  be  characterised  through  the  Rayleigh  quotient: 


If  we  let  ^  **  a  •  Vf#,  w#  can  characterise  the  minimum  positive  eigenvalue 
of  P~lQ: 

tt  aTQa 

x.  =  min 

(m}  oT  Pa 

where  «o  *  (1,0, ... ,0). 

The  following  is  well  known,  see  |l|: 

0  <  **  -  A*  <  C{  inf  i#-xll)’ 

X  —  b  ♦n  for  some  b 
*  €  iVf (A*) 

ii^ii  - 1 


where  A/(A*)  is  the  eigenepace  corresponding  to  A *. 

From  these  observations,  we  conclude  two  things: 

•  The  localisation  of  the  error  estimator  as  defined  in  (4.5)  can  be 
founded  on  exponential  decay  of  the  solution  away  from  "vertical* 
boundaries  and/or  rough  spots  in  the  load. 

•  A  choice  of  basis  functions  is  to  be  preferred  over  another  if  the  first 
leads  to  a  smaller  minimum  positive  eigenvalue  x[*  [nff  *=  0  m  A«), 
sines  such  a  choice  leads  to  the  use  of  less  basis  functions  (a  smaller 
/V<)  away  from  rough  spota).  That  is  evident  from  the  following  ex¬ 
ample. 

There  is  an  orthogonal  matrix  0  such  that  0TP~lQ0  *  D ,  being 
diagonal.  Setting  V  m  0V  yields  the  following  system  of  O.D.E.s 

-V  +  ±DV  -  0T  P~lR  =  C 

with  tht  tolulion: 

M*)  =  x  (vf(*)) 

uf(*)  -  !’ t-V^iG,  +  /*«/“?*(?,</» 

Jo  Jo 
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for  •  >  1.  If  one  for  example  takes  0  »  $(xn),  the  solution  K  involves  terms 
of  sinhfv/^TjJ)  c®*h "here  Hornes  clear  (bat  a 
smaller  x{*  improves  localisation. 

For  N  =  2,  choosing  the  basis  functions  a «  in  (5.1)  yields  xjj  =  0,xj  = 
15,  the  latter  approximating  well  the  eigenvalue  Aj  *  t3  with  respect  to 
exponential  decay.  That  is  also  the  best  one  can  do  given  the  span  for  N 
«  2  (using  (1,  it2}).  In  contrast,  if  one  omits  do  a*  the  first  basis  function, 

B  0  =  *j.  For  N  >  2,  the  approximation  of  Aj  can  not  get  any  worse. 
We  therefore  choose  the  basis  functions  as  in  (5.1). 

Our  initial  computations  suggest  a  practical  confirmation  and  viability 
of  many  of  the  features  described  here  of  this  method.  It  should  be  noted 
that  we  have  not  dealt  with  the  issue  whether  or  not  this  feed  back  method 
is  adaptive,  i.e.  whether  or  not  this  feed  back  method  is  optimal  with 
respect  to  some  performance  measure.  It  will  be  dealt  with  elsewhere. 
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